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Absract—This paper proposes malttimescale voltage 
stability-constrained velWVAR optimization to. improve. the 
Heady-sate stability and security of smart distribution grids, Te 
‘ims to coordinate the pst contingency slow-timescae corrective 
Actions provided by onload tap changer with fastsimescale 
Corrective actions ‘provided. by photovoltaic inverters and. 
battery energy storage system. This multi-timescale operation 
exploits the difference between short- and long-term emergency 
ratings of line flows and voltage magnitudes usually defined by 
utilities for contingency cases. Global optimum Is achieved by 
formulating the peoblem as mixed-integer second-order cone 
program (MISOCP). To justify the investment in battery storage 
the proper coordination between base case slateof-charge and 
Storage dynamic reserve margin required for port-contingeney 
Corrective actions is performed. Also, the steady-state voltage 
Stability margin is computed through detailed voltage stability 
analysis, Simulation results demonstrate the effectiveness of 
proposed multi-imeseale corrective actions in maintaining fhe 
Steady-state stability and security under contingency cases 


Index Terms—Battery storage dynamic reserve margin, convex 
relaxation, distribution grids, multi-timescale corrective actions, 
voltage stabily-constrained volUVAR optimization. 


1 INTRODUCTION 


OLTAGE stabilty-constrained vol/VAR optimization 

(VSC-VVO} model ensures pre- and post-contingeney 
steady-state voltage stability subject to security (eg, line 
fow and voltage magnitude limits) and economie constraints 
[1]. The VSC-VVO model solves the traditional volUVAR 
optimization (VVO) problem while accounting for voltage 
stability limits in base and contingency cases. The growing 
load demand, increasing distributed energy resources (DERS), 
and competitive energy market usually lead to the operation 
of modem distribution grids closer to their operational limits. 
The robustness of the VVO solution against voltage instability 
can only be ensured by considering stability constraints in 
stressed distribution grids (2), [3]. 

Banery energy storage system (BESS) can be used for 
various applications such as energy arbitrage, peak shaving, 
and providing the stabilization required for the integration of 
intermittent renewable sources [4]. The heavy investments in 
BESS for post-contingeney corrective control actions might be 
justified only when the BESS is simultaneously used for oter 
Standard applications during base case, Several researches 
have exploited the fast response capability of BESS [5]-110] 
Reference [S] reduced the required fast-ramping expensive 
generation through BESS. In [6], the photovoltaic (PVbatery 
system is introduced into a security-constzained unit commit 
ment for transmission congestion management. Reference [7] 


employed the system generators and microgrids (inel, storage} 
to relieve postcontingeney overflows, Only microgrids were 
used for corrective control and these were not coordinated with 
generation redispatch to provide multi-timescale corrective ac- 
tions. By contrast, reference [8] implemented multi-timescale 
corrective conirol scheme to relieve post-contingency over 
loads. The bulk and distributed energy storage exploited the 
short-term emergency (STE) and long-term emergency (LTE) 
ratings of lines, respectively. Similarly, the BESS reduced the 
line flows below STE rating and then generation redispatch 
kept the line flows below LTE rating afier contingency [9]. The 
BESS state-f-charge (SOC), which significantly influences 
the storage ability to provide post-contingency corrective ac- 
tions, was not considered, Reference [10] bridged this gap by 
assuming pre-specified initial BESS state-of-charge for multi- 
limescale corrective actions. However, the proper coordination 
mechanism between base case state-f-charge (very important 
for standard applications like arbitrage) and the BESS reserve 
margin required for post-contingency corrective actions was 
not considered. Above literature [6]-(10) mainly considers the 
transmission congestion management task and ignores voltage 
instability issue, which is vital in stressed power grids. 
Voltage stability-constrained voly VAR optimization is a 
mixed-integer nonlinear programming (MINLP) problem. The 
nonconvex equality constraints and presence of discrete con- 
trol variables make it an NP-hard problem. Voltage stability 
margin (VSM) may be indirectly approximated by several 
heuristic voltage stability indices that provide suboptimal 
solutions for computational gains. By contrast, optimization 
approach involving two sets of constraints (one for normal 
and other for critical condition) provides accurate and direct 
means to compute VSM [II]. Reference [12] proposed a 
statistical approximation procedure to calculate VSM. In [13] 
the VSM is directly computed using optimization approach 
for reactive power planning problem. Reference [14] reduced 
the index of proximity to steady-state voltage collapse using 
non-dominated sorting genetic algorithm (NSGA-II). Simi- 
Tarly. particle swarm optimization (PSO) improved the voltage 
stability levels of buses in distribution grids (15]. Although 
evolutionary algorithms (e.g, PSO and NSGA-II) are flexi- 
ble, they may not find a global optimum and usually lack 
performance guarantees due to nonconvex VSC-VVO mode, 
This model may be convexited and solved efficiently through 
convex relaxations. For radial networks, the convex relaxations 
are always exact and yield global optimum provided there 
are no upper bounds oa loads (sufficient but not necessary 
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condition) [16]. [17]. Reference [18] achieved global optimum 
of securty-constained optimal power flow problem through 
semidefinite program. However, it only considered line capac 
ity constraints and ignored voltage stability constraints. Refer- 
ence [2] formulated second-order cone program to reduce gen- 
eration cost in voltage stability-constrained optimal power fow 
(VSC-OPF) problem. Here, the voltage stability constraint was 
approximated by a second-order conic representable voltage 
stability index. However, the BESS was not incorporated and 
multi-timescale corrective actions were not considered. 

The recent transition from passive to active distribution 
networks demands an evolution from offline to online opt 
mization. Reference [19] formulated a linearized power system 
model to solve multistage VVO problem in (near) real 
In (20), the power quality is improved by real-time voltage 
regulation scheme involving smart inverters. Reference [21] 
solved an online optimization problem to ensure the feasible 
operation of distribution grids under forecasting uncertainties 
by proposing a reconciliation algorithm. Our work in [22] 
coordinated the operation of VVO equipment in the presence 
of uncertainties involved in the forecasting of PV output and 
load demand using stochastic programming. Note that the 
forecasting uncertainties have an impact on VVO performance, 
however, the effects of uncertainties on distribution grid op 
eration are beyond the scope of this paper. A comprehensive 
survey of real-time optimal power flow ean be found in [23] 

Traditional VVO equipment such as on-load tap changers 
(OLTCS) and switched capacitors can wansition relatively 
slow (in seconds to minutes) due to their mechanical nature 
while recent DERS” converters can transition very fast (in 
milliseconds) [22]. This situation essentially enables the muli- 
timescale post-contingency corrective actions as the OLTCs 
with slow response and DERS” converters with fast response 
are effective control devices to achieve steady-state voltage 
stability [24]. Huge literature has considered the coordinated 
control of OLTCs with battery storage in standard VVO model 
1254-127]. Reference [25] reduced losses by the centralized 
control of OLTCs and battery storage. Similarly, the coor 
dinated control of OLTCs and storage reduced losses along 
with voltage profile improvement [26], Reference [27] formu- 
lated mixed-integer second-order cone program (MISOCP) 10 
minimize energy purchase cost by coordinating OLTCs with 
battery storage, Only fewer applications of multisimeseale 
VVO models coordinating traditional VVO equipment with 
battery storage can be found in the literature. In our previ- 
ous work [22], we implemented multi-imescale coordinated 
control of OLTCs and battery storage to reduce losses and 
energy purchase cost under forecasting uncertainties. Slow and 
fast timescales were coordinated using two-stage stochastic 
programming. Here, we now consider VSC-VVO model incor- 
porating multi-ümescale post-contingency corrective actions, 
In [28], intrarhour voltage regulation was performed through 
storage by assuming separation in slow and fast timescales. 
However, stability constraints and contingency cases were not 
considered. Reference [29] reduced losses through capacitors 
and inverters operating on slow and fast timescales, Again, 
stability constraints and battery storage were ignored and 
muki-amescale corrective actions were not implemented. 


The motivation behind this work is to suggest a novel multi- 
timescale VSC-VVO model that not only exploits the STE and 
LTE ratings of lines (similar to [8], [9], and [10]) but also takes 
advantage of different steady-state voltage magnitude ratings 
usually defined by utilities for base and contingency cases. 
It is very clear that existing literature has not considered the 
multi-imescale corrective actions provided by the coordinated 
control of OLTC and battery storage. This paper considers the 
post-contingeney fast-timescale corrective actions provided by 
battery storage and PV inverters while OLTC provides post- 
contingency slow-timescale corrective actions. In contrast to 
linearized or DC power flow model in [5], [7], and [19] 
the focus of this paper is to achieve the global optimum 
of original MINLP model by formulating the problem as 
MISOCP. The proposed method offers several technical and 
economie benefits to distribution system operator and cus- 
tomers. These benefits include enhanced system reliability, 
improved system voltage profile, reduced line losses, and 
improved environment due to more DERs' integration. To 
justify the heavy investment in BESS, this paper introduces 
A proper coordination mechanism between base case state-of- 
charge (very important for standard applications like arbitrage) 
and the BESS reserve margin required for post-contingency 
corrective actions. Also, the detailed voltage stability analysis 
is performed to compute voltage stability margin. 

IL, MULTI-TIMESCALE CORRECTIVE CONTROL 

This section presents a conceptual illustration of multi- 
timescale corrective control scheme through simulation. I also 
formulates a preventive corrective VSC-VVO model 


A. Multi-timescale Corrective Control Framework 

Utilities usually not only define various thermal limits (e.g. 
normal, STE, and LTE ratings) to avoid damage to power 
components but also specify several steady-state voltage levels 
EirGrid system maintains 0.95-1.1 pu in base case and 0.9-1.12 
Pu in N-1 contingency cases [30]. Similarly, Westem Power 
Australia allows +6% during normal conditions and +10% 
during emergency conditions in distribution grids [31]. Fig 
1 shows a three-node radial system for simulation purpose. 
Fig. 2 provides the explanation of time axis used in this 
paper. The contingency occurs at instant t; and fast-timeseale 
corrective actions provided by BESS are applied at instant tz 
The ramping of slow-timescale corrective actions provided by 
OLTC is initiated at instant funy and concluded at instant fs 
The STE rating is considered [rom t to ts while LTE rating 
(usually same as normal rating) is considered from fy onwards 

The effectiveness of multitimescale corrective actions in 
mitigating line overflow and vollage limit violations is demon- 
strated by simulation. Fig. 3(a) shows that both line flow and 
voltage exceed STE limits when PV inverter contingency is 
realized at instant ty. The increase in line flow Fia due to 
constant power load results in greater voltage drop across 
line and transformer impedances. Fig. 3(6) shows that BESS 
active power Pzgs increases at instant ta to bring Fy within 
limit and reduces voltage drop. Also, the BESS reactive power 
Quess provides voltage support to bring Vz back to voltage 
STE rating, The tap movement concludes at instant ts and 
both line dow and voltage are within LTE limits. 
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Fig. 3. (a) Voltage and line ow) Power injections and tap poston 
B. Multi-timescale Preventive-corrective VSC-VVO Model 
This subsection integrates the voltage stability constraints 
into the formulation given in [9]. This paper computes voltage 
stability margin through optimization approach involving two 
sets of constraints (one for normal condition and other for 
critical condition). Here, steady-state VSM is quantified in 
terms of MW distance from the operating point to critical 
point or point of voltage instability (PVI) [3] 
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contingency eases and K 
(positive integer) contingencies, The Cia and Cuz denote 
base case scenarios at normal and critical conditions, respec 


tively. For contingency number k € K, the Cy and Ck, rep 
resent the kth fast-imeseale contingency cases at normal and 
vical states, respectively. Similarly, the Ath slow-imeseale 
onungency cases at ommal and critical conditons are de. 
seribed by Cy and Chp, respecively. The multitimescale 
preventive-cortective VSC-VVO model can be formulated as 
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control variables (eg, BESS activafcactive power and PV 
inventer reactive power output) anu" defines slow-timeseale 
contrals e'g.. OLTCS tap positions) for scenario e: Tae 2° 
describes siale variables and o°°) refers to probability of pre 
and post-coningeney scenarios, The equality constraint set 
10a). 38), (48), (Sa), 69), ail. mainly describes Joad ow 
equations. Equation 2a) also includes BESS sae-f-charge 
“The inequality constraint set (2b), (4b), (6b)] sacs phys- 
«al limits of components and operational limits of distribution 
grids at normal condition. In contrast, the inequality con- 
Saint set (3b), (Sb), (7)] mainly includes vllage stability 
constraints involving Joad multiplier at critical condition. 
Depending on the context, the Annez denotes the normal, STE, 
or LTE rating of power component. Note that STE rating is 
considered fram t to ty so constraints associated with ramping 
period are unnecessary and post-contingency slow-timescale 
Constraints are defined only at instant ts. The applicability of 
these constraints and variables are summarized in Fig, 2. De 
tiled explanation af constraints is included in the next section. 

Coupling constraint (8) mits the movements of contral 
variables between base and_post-contingeney cases. This 
paper assumes unrestricted variations of control variables 
(Auf mas =A ages) IF post-contingency constraints 
ane violated then VSC-VVO model firs takes the corrective 
measures to ensure security and stability. This is due to 
Small cost of corrective actions as the contingencies occur 
taely and do not last for Tong time period: The preventive 
measures are taken only if necessary as the preventive actions 
{a\2*? aS") may operate the system away from optimal 
posal for long period respective ol contingency event This 
makes preventive actions costlier than comeenive actions [32] 

Ul. PROBLEM FORMULATION 

This section converts the original MINLP problem to MIS- 
OCP model to achieve global optimum. Proper coordination 
mechanism for base case state-of-charge and BESS dynamic 
reserve margin required for corrective actions is presented 
A. MINLP Formulation for Mal-rmescale VSC-VVO Model 

Fig, 4 shows a radial distribution network with set of buses 
N= {0,1,2 Nuus} and its oad flow equations can 
be represented by DisFlow branch equations (38). Let pe 
and quy be et active and reactive power injections at bus 1 
respectively. Define Piya, Qa amd Lay as the active power, 
teacdve power, and squared current flows ftom bus i to f, re- 
spectively. The voltage magnitude and squared voltage at bus $ 
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ig 4. Single line dara of radial aon system [22.133 
are denoted by s,s and Uy, respectively. The objective of pro- 
posed muli-timescale VSC-VVO model is to minimize losses 
(9) while ensuring voltage stability and security constraints. 
The mixed-integer nonlinear se model is given by 
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for every CC and t © 7, where Tie test of time instants 
upto Ty. The ru; and zuy denote resistance and reactance of 
branch (i, j), respeeiveiy. Tae tapy is ap position of OLTC, 
Nap refers to total number of taps, and 927 describes ca 
pacitor reactive power The PV active and reactive powers are 
denoted by p, and g7, while active and reactive load demands 
are indicated by pf, and of, The BESS charging, discharging, 
and reactive powers are denoted by p74, pPI, and q2 553, 
respectively. The apparent power rang of PV inverter and 
line (3) are indicated by 8,0 mar 888 Sij lineer 


Bas. (10}+12) describe DistFlow branch equations for lines 
while (13) indicates DistFlow equations far substation OLTC. 
Eqs. (15)-(16) assume plusminus signs for power injec- 
Uon/extrction, When certain power equipment is not installed 
at a bus, corresponding entry will become zero in (15}-(16). 
Section II-B includes detailed BESS model, The limits on 
taps, voltages, and apparent powers are indicated by (1720). 


Multi-timeseale corrective actions are realized by (21}(22) 
The contral setpoints at critical condition are assumed to be 
same as setpoints at normal condition for computing voltage 
stability margin, The slow-timescale controls provided by 
OLTCS remain same between base and fast-imescale contin- 
agency cases. This exhibits preventive behavior (21a) during 


ppost-contingency short-term period and corrective behavior 
(21b) during post-contingency long-term period. By contrast, 
‘equation (22b) shows that fast-imescale corrective actions 
provided by PV units (also by battery storage) are applied 
during post-contingency short-term period 
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Consider © € (Cpu. Chr: C$n] and k € K. Let pf, be load 
power factor and V'SMiarset be the pre-specified VSM target, 
The difference in the “normal” and “eritical” states’ total MW 
load is kept apart by a minimum VSMorgey [12]. Note that 
voltage stability margin constraints (23) [13] are included only 
al critical states. This paper assumes proportional increase of 
loads (23a) and maintains constant load power factors (235) 
st critical states, Other loading directions ean be considered 


‘based on the expected pattern of increase fron the initial loads 
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IEEE task force on val/VAR contro! mentioned PV inverter 
shutdown event [35]. This paper assumes all N-1 PV inverter 
outage contingencies. Line outages are not considered because 
any fault in the grid could result in islanding and disconnection 
‘of loads. The outage of PV unit at bus i can be modeled as 
ag D oo 
‘The MINLP model (9-24) for miltsimescale VSC-VVO is 
nonconvex due to quadratic equalities in (13b) and (14) 


B. Battery Energy Storage System 

‘The batery storage model is given by (25) [17] with c € C 
The dynamic equality constraint (25a) is convex and links ad- 
jacent intervals of duration At. Equal initial and final state-of- 
charge are ensured by (25c). Eqs. (25b)-(25e) are linear con- 
straints while (259) is convex quadratic inequality constraint. 


SOCS = o SOCEM 

4 AE (m0HA OHNE) pIE PIS) sy) 
04 SOCEM < SOC mae as) 
socn ~ soog (250) 
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where SOC, is the storage state-of-charge at instant t and 
SOC, mas iS energy capacity rating. The on, 18", and nP7S 
denote self-discharge. charging, and discharging efficiencies, 
respectively. The BESS charging, discharging, and apparent 
power ratings are respectively defined by PCH Amas 
PDISumar: and Sinessmar: The binary variable fos 
prevents simultaneous charging and discharging of storage. 
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Define C* = (Cy, Cr: Chw, Ch) as the set of th 
posteontingency scenario, where C% © C, c € C4, and 
K € K. This paper introduces a proper coordination mech- 
anism between base case state-of-charge (very important for 
standard applications like arbitrage) and the BESS reserve 
margin required for postcontingency corective actions by 
replacing (25h) with (26), The BESS reserve margin includes 
the amount of energy that must be delivered or stored to guard 
against all the contingencies. I also includes the BESS cont 
‘bution towards the reduction of postcontingeney cost function 
weighted by contingency probabilities. Note that maximum of 
conver function is convex so it can be seamlessly integrated 
ino inequality constrain without compromising its convey. 
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26) 
Previous work sich as [9] suggested similar reserve margin 
constraints (w/o 1/74 and P15) but didn't coordinate them 
with BESS state-of-charge. By contrast, our dynamic approach 
is less conservative and more economical. Note that slorage 
systems are usally seen as active power sources and addi. 
tonal reactive power capability comes from power electronics 
ase to connect them to distribution grids. The constraints (19) 
snd (258) governing reactive power capability of PV units and 
Storage are sutheient for the purpose of VSC-VVO model [17] 
E MISOCP Formulation of Matt-timercale VSC-VVO Model 
This subsection presents MISOCP formulation for multi- 
timescale VSC-VVO model. Reference [33] used binary ex 
pansion scheme to replace (17) by (27a) and big-M method to 
transform nonconvex constraint (13b) into mixed-integer lincar 
expressions, Consider © € C and substituting new variables 
mbar U, 201-20, UG), and yI, ma in 
(1b) and (13) would resul in (270) (27g) ia big M method. 
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where Ane Is a binary variable, Nory denotes binary 
expansion length. My and Ma are big-M method constants 
Exact linearized OLTC model (27) 133] bas replaced non 
convex constraint (13b). Reference [16] describes two relax 
‘on steps to convexify (14). Angle relaxation. step results 
in (12) and (14) by taking squared complex magnitudes, 
Applying second-order cone relaxation step on (14) gives 
L 2 (RE) + OSEAS a 


subject to (10)-(13a), (130) 
‘The global optimum would attain the equality of (28) if 
MISOCP relaxation is exact. For radial networks, the con- 
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vex relaxations are always exact and yield global optimum 
provided there are no upper bounds on loads (sufficient but 
not necessary condition) [16], [17]. The proofs relevant to 
the exactness of convex relaxation for radial grids are given 
by ( [16], Theorems 1 and 4) and ( [17], Theorem 3.1). 
Exact relaxation results in feasible (hence globally optimal) 
solution while inexact relaxation yields infeasible solution that 
only gives a lower bound on true optima [17]. Note that the 
simulation studies in Section V-A and V-B always resulted in 
alobal optimum by attaining equality of (28) 
IV. MULTI-TIMBSCALE MPC FRAMEWORK 

Fig. 5 shows the proposed mulictimescale operational 
framework based on model predictive control (MPC) alge- 
rithm. The advanced metering infrastructure (AMI) system re- 
ceives the equipment measurements and status while sends the 
contral signals to battery storage and PV inverters. Similarly, 
the substation remote terminal unit (RTU) moniter and control 
the tap position of OLTC. The proposed multi-timeseale MPC 
‘framework is dictated by the distribution management system 
(DMS), which is consistent with [19]. Here, the DMS software 
inventory includes distribution SCADA, PV production fore- 
cast tol, demand-forecasting tool, VSC-VVO optimization 
engine, and remedial action system (RAS). Also, the details 
about the AMEDMS integration can be found in [36] 

‘The demand and PV output profiles are assumed to be 
forecasted based on historical dala and actual measurements 
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at instant ¢ 1. These forecasted profiles are used by VSC: 
VVO optimization engine to compute the base case and 
post-contingency corrective set-points for the entire prediction 
horizon between £ and t+ T. The proposed MPC framework 
only implements the frst decision variables and the same 
procedure is repeated for the next prediction horizon. Note that 
the proposed control framework can be seamlessly integrated 
into the day-abead and hour-ahead energy markets. 

The localized contol of PV inverters and OLTC equipped 
‘with line drop compensator (LDC) could result in hunting phe- 
nomenon among equipment that may interfere with network 
coniro. Also, the localized control methods cause increased 
tap operation of OLTC under high PV penetration [25]. By 
contrast, the proposed MPC framework based on centralized 
approach avoids hunting phenomenon and mechanical wear 
and tear due to the frequent switching of OLTC. The proposed 
centralized implementation would be relatively fast as 3t does 
not introduce intentional time delay between 30 seconds and 5 
minutes to prevent hunting among equipment. The tap changer 
equipped with LDC experiences unintended automatic raise 
‘operation due to long-term voltage instability issue [24]. The 
proposed approach enables the OLTC blocking operation to 
avoid the negative actions of OLTC in critical conditions. 

V. SIMULATION RESULTS AND ANALYSIS 

The proposed multi-timescale voltage stability-constrained 
VolUVAR optimization is examined on modified IEEE 33-bus 
feeder [34] as shown in Fig. 6. The BESS and PV units provide 
Past-contingeney fast-timescale corrective actions while OLTC 
provides post-contingency slow-timescale corrective actions. 
Assume that substation voltage vwy is 1.0 pu and OLTC has 

10% tap range and 32 tap positions (ie. tap=[-16,16). 
Previous literature such as [8}-[10} has considered only single 
operating point. This paper assumes PV generation profile 
fom NREL [37] and hourly load profile from SCE [38] 
fon 25th April, 2018 as the computation of dynamic reserve 
margin of battery storage demands mul-period optimization. 
The loading factor in Fig. 7 is the percentage value of init 
‘operating point and highlighted portion is the period of inter- 
est. AlI N-L PY inverter outage contingencies are considered 
and corresponding probabilities are indicated in Table L Here, 
mean time to failure of each PV inverter is assumed to be 800, 
(see [39] for detailed probability calculations), GAMS/CPLEX 
[40] is used to solve MISOCP problem on ordinary PC. 

A VSCWO Model Operating m Fully Corrective Mode 

‘This subsection explores the situation when multi-timescale 
VSC-VVO model operates in fully corrective mode. Initial 
operating point load is 3.715 MW and 23 MVAR. Assume 
that each shunt capacitor has 30 KVAR rating. The BESS 
installed at bus 33 is assumed to have a capacity of 1.1 MVA/S 
MWh, Also, the PV inverters each with the rating of $50 KVA 
are installed at bus 11, 12, 22, and 29. The proposed model 
maintains 0.95-1.05 pu in normal condition and 0.93-1.07 pu 
in STE condition. The STE branch flaw limit (1.56 MVA) is 
assumed to be 1.2 times larger than normal limit (1.3 MVA). 
The LTE ratings for both voltage and line flow are assumed 
to be equal to normal ratings. The voltage 
line flow limits are unbounded at critical c 
operational range can be considered based on system's data. 
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CONTINGENCY SCENARIOS AND PROBABILITIES 
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Fig. 8(@) shows that tap position is high between 14:00- 
15:00 (case Cry) due to increased load demand and reduced 
solar generation. Fig. 8) shows that the battery storage is 
charging between 12:00-13:00 (case Cay) when load demand 
is low and PV output is high. By contrast, the BESS reduces 
power losses by discharging at instant 15:00 (case Can) 
When demand is highest and PV output is lowest. Fig. 8c) 
depicts the voltage support provided by BESS reactive power 
tions for all pre- and postcontingency cases. Note that 
imescale corrective actions are (almost) same for con- 
lingeacies C! and C? due to PVII and PV12 being installed 
ät adjacent buses. The contingencies C? and C denote PV22 
and PV29 outages. Fig. B(d) shows that BESS state-of-charge 
meets the constraint (25e) of same initial and final values 
Fig 9(a) shows the substation active and reactive power fow 
for base case. Active power ow is maximum al instant 15:00 
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due to high load demand and low PV generation. Fig. 96) 
depicts the dependence of capacitors’ reactive power injections 
on voltage magnitudes at the respective buses. Fig. 9(¢) shows 
the total losses including storage charging and discharging 
losses (n°#—y!9=0.05). Reactive power support provided. 
by PV units to guard against all contingencies is shown in Fig 
10. Different reactive power output indicates the reliance of 
corrective actions on PV units’ locations 

Fig., 11 includes the solutions provided by optimization 
model (for normal, short-term, and long-term periods) and 
OpenDS load flow package (for transitional period). Here, we 
assume the time-axis Ge., (a, fa,ta) definitions as described 
in Section IL Fig. 11) shows the violation of line fow 
STE limit during transitional period at instant 15:00 (f). 
The corrective actions then keep the post-contingeney line 
flows within limits. Similarly, the voltage decreases upon the 
occurrence of contingency at instant t1 in Fig. 11b). The 
multi-timescale corrective actions then improve the minimum 
Voltage profile to reduce post-contingency losses 


Fig. 12 demonstrates the proper coordination between base 
case stale-of-charge and BESS dynamic reserve margin re- 
quired for post-conlingency corrective actions. These correc- 
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tve actions contribute towards the mitigation of constraint vie- 
lations and reduction of cost function (e.g., losses). The proper 
coordination mechanisms is provided by constraint (26), Here, 
the corrective actions are provided in the form of discharges 
with minimum overall reserve margin for contingency C and 
maximum overall reserve margin for contingency C“ 

Fig. 13 shows the P-V curves computed for incremental 
Toad increase with fixed control variable settings using quasi- 
dynamic simulations in OpenDSS. The P-V curves are plotted 
al instant 15:00 and bus 18 is assumed to be a representative 
bus due to minimum voltage magnitude, Any other represen- 
tative bus can be considered. The 20% VSM target requires 
Joad multiplier 2 to be at least 125 pu as per (23c). Note that 
operating point at eritical state can be anywhere between mini- 
mum 9 value and point of voltage instability. Here, the critical 
operating point is maintained at 1.25 pu due to the nature 
of post-contingeney cost function (Le. loss reduction). The 
“VSM is reduced from 3.83 pu (case Cau) to 355 pu (instant, 
t1) upon the introduction of C contingency. Muli-timeseale 
corrective actions then increase the VSM to 3.865 pu. Only 
corrective measures were taken to ensure security and stability 
as both VVO and VSC-VVO models resulted in exactly same 
overall loss of 66.02 kW. Aso, the load fow solution failed to 
‘converge beyond nose point and operating points above PVI 
usually represent the satisfactory operating region [3]. 

B. VSC-VVO Model Operating in Preventve-Correctve Mode 

This subsection explores the situation when multi-timescale 
VSC-VVO model operates in preventive-corrective mode. As- 
sume initial operating point to be 22.29 MW and -13.8 MVAR 
with six times increase in constant leading power faclor loads. 
The BESS with capacity of 3.3 MVA/IO MWh, shunt capac- 
itors each with 600 KVAR rating, and PV inverters each with 
the rating of 1.1 MVA are considered. The tap positions are as- 
sumed between 0 and 47, The proposed model maintains 0.92 
1.08 pu in normal condition and 0.9-1.1 pu ia STE condition 

The VSC-VVO model operation in preventive-coreetive 
mode results in solutions different from standard VVO model. 
This is due to additional costly preventive measures taken to 
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meet the security and stability constraints. Fig. 14(@) shows tap 
positions and Fig. 14) reveals the rising tend of BESS active 
power injections (case VVO and Cay) with steep rise at instant 
15:00 due to high demand and low PV generation. By contrast, 
the BESS active power becomes zero at instant 15:00 (case C*) 
owing to very small post-contingency probabilities. Fig. 14(c) 
depicts BESS reactive power flows. The BESS with 50% initial 
state-of-charge is completely discharged to reduce losses and 
constraint (25) is not imposed here as indicated by Fig. 14) 


‘The preventive measures could move the system away from 
optimal operating point as evident from Figs. 15 and 16, The 
substation active and reactive power flows are shown in Figs 
15(a) and 15(b). The capacitor reactive power in Fig. 1$) is 
highest at instant 15:00 (case cap32-Can) due to high voltage 
magnitude at bus 32. Fig. 15(d) shows the higher power loss 
in base case Caw as compared to VVO solution. 


Fender voltage poi at instant 15:00. 


‘The multi-timeseale corrective actions exploiting the short- 
term emergency ratings are shown in Figs. 16 and 17. The 
line flows are depicted in Fig. 16(2) Fig. 16(b) shows that 
voltage LTE limit (ease VVO) is violated between 15:00-16:00 
due to absence of corrective actions in standard VVO model 
By contrast, no such violation is observed (ease VSC-VVO} 
because of coordination beween fast- and slow-timescale cor- 
rective actions. Fig. 16c) shows that multi-imescale approach 
exploits the lower STE limit at instant 15:00 (from tz to t) 
Fig. 17 shows the entire feeder voltage profile at instant 15:00. 
‘The voltages remain within limits and clear rise at profile tail 
(case Cy) is due to high active power injection from slorage. 


ig. 18 demonstrates the coordination between base case 
atatecotcharge and BESS dynamic reserve margin required 
for post-contingeney corrective actions as per (26). The reserve 
margin is nil at instant 15:0 (case C") because the storage dis- 
charges completely to reduce base case loses, Fig. 19 shows 
the P-V curves with VSM{CE4) =49% and VMI, =10% 
‘The voltage stability margin at optimal VVO solution i 40% 
so VSC-VVO model moves the voltage upward from 0967 pu 
10 0.9985 pu to meet 43% target This preventive operation re- 
sults in suboptimal solution. The overall osses are 5 0147 MW- 
VVO) and 5:325 MW (ease Cox) indicating 6.18% increase, 


The proposed MISOCP model yields global optimum by 
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attaining equality of (28) with zero optimality gap. Figs. 3-19 
clearly indicate the feasibility (hence optimality) of solution 
returned by exact MISOCP relaxation. We have also performed 
a separate optimization to maximize # subject to fixed contol 
setings using GAMS. The computed point of voltage instabil- 
ity by Open DSS and GAMS is found to be same indicating the 
accuracy of results, Simulation times for VVO and VSC-VVO 
in fully corrective mode are 1.83 s and 444.43 s, respectively. 
VI, CONCLUSION 

‘This paper coordinated the post-contingency slow-timescale 
corrective actions provided by OLTC with fast-timeseale cor 
rective actions provided by BESS and PV units. This muli 
timescale operation exploited the dilerence between STE and 
LTE ratings of line Rows and voltage magnitudes. Unlike the 
standard VVO model, the proposed VSC-VVO model can 
maintain the voltage stability and security under contingencies, 
MISOCP model reduced losses and achieved global optimum, 
Also, the preventive actions are shown to move the stressed 
system away from global optimum. Proper coordination be: 
tween base case stateof-charge and storage reserve margin 
required for post-contingency corrective ations is ensured. 
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